CERN-TH.7164/94 
SWAT/93-94/26 
|hep-th/9402084| 
February 1994 



The exact mass-gaps of the principal chiral models 



Timothy J. Hollowood* 

CERN-TH, CH-1211 Geneva 23, Switzerland. 
hollow@suryal 1 . cern . ch 



ABSTRACT 

An exact expression for the mass-gap, the ratio of the physical particle mass to the 
A-parameter, is found for the principal chiral sigma models associated to all the classical 
Lie algebras. The calculation is based on a comparison of the free-energy in the presence 
of a source coupling to a conserved charge of the theory computed in two ways: via the 
thermodynamic Bethe Ansatz from the exact scattering matrix and directly in perturbation 
theory. The calculation provides a non-trivial test of the form of the exact scattering 
matrix. 
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1. Introduction 



The principal chiral, or QxQ, sigma models are two-dimensional quantum field theories 
that are integrable at the quantum level. The fact that the theories are integrable means 
that their scattering matrices are factorizable. Such S-matrices have been conjectured 
for all the theories corresponding to the classical Lie algebras [1]. Expressions for the 
complete S-matrices for SU(iV) can be found in [1] and for Sp(2iV) in [2]. Not all the 
S-matrix elements are known in the case of SO (AT). It is an outstanding problem to prove 
from first principles that the S-matrices actually do describe the lagrangian field theories. 
This is especially important because a given factorizable S-matrix is ambiguous since it 
may always by multiplied by CDD factors. Connecting the S'-matrix picture with the 
lagrangian picture is highly non-trivial since the masses are generated dynamically and 
the theories are asymptotically free. 

In a series of papers such non-trivial tests have been applied to various integrable 
models: the Gross-Neveu model [3], the O(N) sigma model [4,5] and the SU(iV) principal 
chiral sigma model [6,7] (the SU(2) case was also considered in [8,9]), using a technique 
known as the Thermodynamic Bethe Ansatz (TBA). The central idea is to couple the 
theory to a particular conserved current and then compute the response of the free-energy 
for large values of the source in the regime when conventional perturbation theory is valid. 
The same quantity can then be computed directly from the S'-matrix using the TBA 
equations at zero temperature in which the coupling to the source appears as a chemical 
potential. By comparing the two expressions a non-trivial test of the S-matrix is obtained 
as well as an exact expression for the mass-gap (the ratio of the physical mass to the 
A-parameter). 

In general the solution of the TBA equations at zero temperature coupled to an 
arbitrary chemical potential would be a formidable problem, even in ultra-violet limit, 
since the equations are a set of coupled integral equations. However, by a judicious choice 
of the source the state of the system contains just one particle which undergoes elastic 
scattering (the particle is the highest weight state of a multiplet). The TBA equations 
then reduce to a single integral equation which can be solved in the ultra-violet limit using 
generalized Wiener-Hopf techniques [4,10] (for a summary see the appendix of [3]). 

In this paper we extend the results of [6] to the principal chiral models for all the 
classical Lie algebras and arrive at a universal formula for the exact mass-gap. We also 
show that by tuning the source we can force the system into inequivalent ground-states 
which each consist of a single type of particle. The fact that the ground-states are pure 
for particular values of the source, is presented as a conjecture whose ultimate justification 
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comes from the agreement with perturbation theory; however, it should be possible to 
prove this fact directly from the full TBA equations of the models. 
The principal chiral models are described by a lagrangian density 

C = -^TT(g- 1 d^g.g- 1 d^g), (1.1) 

where g is a group valued field. The theory is invariant under a global symmetry cor- 
responding to left and right multiplication by the group g i— > h^gh^ , /il,r £ Q- A is a 
dimensionless coupling constant. 

The S'-matrices that have been conjectured to describe the scattering of the states 
of the model describe r = rank(^) particles which are associated to the fundamental 
representations of Q. The masses of the particles, except for those associated to the spinors 
of SO(iV), can be described by the universal formula [1] 

m a = m .\ ;7 , a = 1,2,..., 1.2 
sin(7r/#) 

where g is the dual Coxeter number of the Lie algebra associated to the group: for A r , B r , 
C r and D r it is r + 1, 2r — 1, 2(r + 1) and 2(r — 1), respectively. The masses of the spinors 
of B r and D r are 

777 777 

B r : m r = D r : m r -i = m r = (1.3) 

2sin(7r/fif) 2sin(7r/fif) 

The particle with mass m a transforms in the following representation of Q x Q [1] : 
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(1.4) 



where V a is the a th fundamental representation of Q with the standard labelling of the 
Dynkin diagram [1]. Notice that although the particles are associated to the fundamental 
representations they are sometimes reducible in the case of SO(iV). 

Fortunately, we shall not require the expression for the complete S'-matrices but only 
those elements for the particles of the highest weight in each multiplet (so with quantum 
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numbers |u; ,u; a ) where the w a 's are the fundamental weights). The S'-matrix amongst 
these states is purely elastic and their expressions can be extracted from [1] : 

Sab(O) = exp liirSab + 2i J — sm(0x)[R ab (x)-6 ab ]j, (1.5) 

where 9 is the rapidity difference of the incoming particles and the kernel R a b(9) has the 
following form for all the particles except the spinors: 

n • i / min(a,6) \ • u / maxfa.M \ 

2 sinh — ^r^nx sinh . v — nx 



sinh (^ra) 



B r , C r , D r : Rab(x) = 



(1.6) 

2 sinh I ^ 2 - J -i:x j cosh I 2g — ttx I 

COsh (^7Tx) 



The S'-matrix elements involving the spinors can also be deduced from the formulas of [1] 
but we shall not require them. 

In the following two section we calculate the free-energy in the presence of a source 
coupling to the conserved charge of the Q x Q symmetry in two ways: from the lagrangian 
using perturbation theory and from the ^-matrix using the thermodynamic Bethe Ansatz. 



2. Free-energy in perturbation theory 



The conserved currents of the left and right symmetry are = g~ 1 d^g and = 
(d fx g)g~ . We wish to couple to modify the hamiltonian of the theory by introducing a 
coupling to the conserved charge of the diagonal action of the symmetry. At the lagrangian 
level this is described by introducing the "covariant derivative" [6]: 

D^g = d^g -ihS lx0 (Qg + gQ), (2.1) 

where Q is a constant element of the Lie algebra. The lagrangian density in the presence 
of the source is 

£ = £o- -^Tr ((g-'Q + Qg~ l )d Q g) - -^Tr (Q 2 + g^QgQ) . (2.2) 

The quantity we will calculate is Sf(h) = f(h) — f(0) where f(h) is the free-energy 
per unit volume in the presence of the source. We shall perform a perturbative calculation 
in the running coupling X(h) which in the ultra-violet regime (large h) runs to zero and 
hence is the regime where perturbation theory will be reliable. We shall only perform the 
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computation to one loop; however this will be sufficient to provide a non-trivial check of 
the S'-matrix and allow for the evaluation of the exact mass-gap. 
An explicit basis for g is provided by 



9 = exp 



jz^n (a) £; Q + m-#J , (2.3) 



where the fields satisfy the reality condition n^ a '* = n*- - "-* and n* = n, and the sum is 
over all the roots of algebra. In the above E a is the usual step generator associated to 
a root a and H is the generator of the Cartan subalgebra. In what follows we choose a 
normalization in which the roots of the simply-laced algebras have length-squared 2 and 
the the long roots of B r and C r have length- squared 2 and 4, respectively. 

Without loss of generality we take Q to be in the Cartan subalgebra so Q = q ■ H, 
where q is some r-dimensional vector. The quadratic part of the (euclidean) lagrangian 
density (2.2) is simply 

Ah 2 



C = 



q 2 + ^J2 {d^O^-^ + h 2 (a- qfn^n^ } , (2.4) 



A 2 

a>0 



where the sum is over the positive roots and for simplicity we have changed the normal- 
ization of some of the n^'s. Notice that the Cartan subalgebra fields n are completely 
decoupled to this order in the loop expansion. 

The tree level contribution to Sf(h) is simply 

Ah 2 

Sf(h)o = --^q 2 . (2.5) 

To evaluate the one-loop contribution we use dimensional regularization. Using standard 
methods one finds 

5f{k)l = "S^ 2 + £ E(« • I 1 " 7E + hi47r - In (h 2 (a ■ q) 2 / > 2 ) } + •••, (2.6) 

a>0 

where e = d — 2, fx is the usual mass parameter of dimensional regularization and g is the 
dual Coxeter number as before. To cancel the divergence in the MS-scheme we add to the 
lagrangian a counter-term 

The quantity 5f(h) is renormalization group invariant when A runs with \x. We can use 
this freedom to set \i = h. The way that the coupling constant runs with h is determined 
from the form of the counter-term. One finds 

h^X 2 = -£\'-(3 1 X 6 -0(X% (2.8) 
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although the second universal coefficient of the beta-function (3\ is not determined at the 
one-loop level. The expression for the first coefficient of the beta- function f3 = g/8n agrees 
with [11]. 

The expression for the free-energy is then 



Ah 2 h 2 

' ' a>0 



where the explicit h dependence is obtained by expressing the running coupling in terms 
of the A-parameter by solving (2.8): 

= (3 In + ^ In In ^- + O { —^—\ , (2.10) 
\ 2 {h) M0 Aj^ Po Ams V ln A^7 

where f3 = g/8n. Equation (2.9) is the generalization to all the classical Lie algebras of 
equation (17) of [6] for A r . 

For comparing with the expression for the free-energy from the TBA calculation we 
set q = u) a /(2uj 2 ) (excluding the spinors of SO(AT)). Writing 



Sf(h) = -^kl 



h A Pi , , h ^ 1 

■ J ' J- ^lnln-r— + - h 

l MS Po A MS \ln^z^ 



(2.11) 



By explicit computation we find for A r that 



K ~ 2na{r + 1 - a) ' Aa ~ [n \2a(r + 1 - a) J 2' ^ L2) 
and for the other algebras a universal form applies: 

kl = -^, A a = -lna-l-^^ln2, (2.13) 
2na 2 g 

where the quantity d\ is the dimension of the vector representation of the algebra, i.e. 
r + 1, 2r + 1, 2r and 2r for A r , B r , C r and D r , respectively. 



3. Free-energy from the .S-matrix 

In this section we will calculate 5f(h) in the ultra-violet limit, h 3> m directly from 
the S'-matrix. The technique is known as the Thermodynamic Bethe Ansatz (TBA) and 
in its most general form it allows one to calculate the behaviour of the free-energy of a 
one-dimensional gas of particles described by a factorizable S'-matrix on the temperature 
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and in the presence of a chemical potential. The free-energy is given in terms of a set of 
functions — in general infinite in number — which satisfy a set of coupled integral equations 
(the TBA equations). 

For our application we working on the plane and hence at zero temperature. The 
coupling of the theory to the source in (2.2) leads to a particular form for the chemical 
potential. The one-particle states are labelled by two weight vectors v) (as well as 
the rapidity) and they can be chosen to be eigenstates of the charge Q with Q\n,i>) = 
q ■ (p + v)\n>-> v). The full TBA equations for the principal chiral models are known [1]; 
however, for particular choices of q, extending the philosophy of [3-6], we conjecture that 
only one particle contributes to the ground-state and the infinite set of TBA equations 
reduces to a single equation. The precise formulation of our conjecture is that when 
q = uj a /(2uja) only the unique particle with the highest charge/mass ratio contributes to 
the ground-state, i.e. the particle \oj a ,ui a ) which is highest weight state of the multiplet 
W a . This particle has Q eigenvalue 1. However, we exclude the the spinor particles of 
SO (AT) from this conjecture. We shall find that this proposal leads to a result which is 
perfectly consistent with the perturbative calculation. Notwithstanding this, it should be 
possible to prove the conjecture directly from the full TBA equations. 

The expression for the free-energy with q = u a /{2uj'^) is then given in terms of a 
quantity e(6>) which satisfies the integral equation: 



J-B 

The parameter B is determined by the boundary condition e(±B) = and the kernel is 
given by 



where S aa (9) is the S'-matrix element of the particle \uj a ,u) a ) with itself (1.5). Once e(6>) 
is known the expression for the free-energy per unit volume is 



Our problem is to solve the integral equation (3.1). In general it is not possible to 
find the solution of such an equation in closed form; however, for comparing with the 
perturbative result we only need to compute the free-energy in the ultra-violet regime h ^> 
to. In this limit a series solution can be found using generalized Wiener-Hopf techniques 
[3,4,10]. The first problem is to decompose the kernel R aa (x): 




(3.1) 




(3.2) 




(3.3) 



Raa(x) 



1 



(3.4) 
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where G±(x) are analytic in the upper/lower half-planes, respectively, and G_(x) = 
G+\—x). The next step in the solution technique depends upon the form of G+\x). For 
all the principal chiral models 



G^m = ^={i-b a ^o{e)}, 



(3.5) 



for constants k' a and b a . So in this respect these models are of similar type to the O(N) 
sigma model rather than the fermion models. With G+(i£) of the form (3.5), [6] gives a 
formula for the first few terms in the expansion of the free-energy 



Sf(h) = 



, h 1 

In h In 

m„. 



G<?>(i) 



l + \\n\n— + [ —Lr- 
2 m a \ In -A 



""MS 



(3.6) 



The explicit expressions for the decompositions are for A r 
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(3.7) 



For the other algebra one finds the universal form 



g 

2yra£ 



r (i + it) 

f . / a . a o-2a, q — 2a 1 1 

x exp ^ -M - In - + In In - 

I V<7 <7 2^7 2^7 2 2 



(3.8) 



From these expressions we find that k' a equals k a in (2.12) and (2.13) and for A r 

, r + 1 — a , r + 1 — a a a 
6 a = — — in — : 1 in 



r + 1 



r + 1 r + 1 r + 1' 



(3.9) 



whilst for the other algebras 



a ln a ^-2a ln ^-2a_ l_ ln l__ 2a 
a g g 2g 2g 2 2 g 



(3.10) 



Comparing the expression (3.6) with the result of the perturbative calculation (2.11) 
we see that they are in perfect agreement if m a oc sin(7ro/fif) which is true for all the 



particles excluding the spinors of SO(iV) (1.2), and furthermore the expression for the 
mass-gap has a universal form: 



9 ex V l( 2 AM)^\^(-), (3-11) 



A MS V 7 ^ (A 2 9 

where m is the mass of the vector particle, which is the lightest particle in the theory (since 
without loss of generality it is only necessary to consider B r for r > 3 and D r for r > 4). 

In addition the S'-matrix calculation implies that the universal ratio Pi/ Pq = 1/2 in 
exact agreement with the perturbative calculation of [11] a fact first pointed out for the 
SU(r + 1) theories in [7,8]. The expression for the mass-gap (3.11) reduces to that of [6] 
for A r . The explicit expressions for each group/algebra are 



SU(r+ 1), A r : 


m 

A MS 
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= ^ 2 

y/iie 
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(3.12) 



2r - 2 



The fact that the perturbative result and the S'-matrix result are consistent provides 
strong grounds for believing that our conjecture about the structure of the ground-states 
is correct. As has been pointed out in [6], the fact that the TBA calculation reproduces 
the universal part of the beta-function Pi/Pq is a highly non-trivial fact. In addition if the 
S'-matrix were modified with CDD factors then the thermodynamics would be drastically 
altered and the perfect agreement with the perturbative result would be destroyed. It 
would be interesting to compare these results with lattice simulations. 

I would like to thank Tim Morris and Michel Bauer for very useful discussions. 
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